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We present a generalized integral fluctuation theorem (GIFT) for general diffusion processes using 
the Feynman-Kac and Cameron-Martin-Girsanov formulas. Existing IFTs can be thought of to be 
its specific cases. We interpret the origin of this theorem in terms of time-reversal of stochastic 
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I. INTRODUCTION 



Feynman-Kac (FK) formula, originally found by Feynman in quantum mechanics [l| and extended by Kac 0], 
^ . establishes an important connection between partial differential equations (PDEs) and stochastic processes. Using 
^ ' this formula, one may solve certain PDE by simulating a stochastic process. Recently, the remarkable FK formula 
. was employed by Hummer and Szabo (HS) Q to prove the celebrated Jarzynski equality ( JE) 0, [HI in nonequilibrium 
I ' statistics [6| . Their work not only provided a concise proof but also pointed out that intrinsic free energy surfaces of 
"j^ , biomolccules could be efficiently extracted by nonequilibrium single-molecule manipulation experiments. 

' Compared to extensive interests of applications of HS's work in single- molecule biophysics 0,13, little attention was 
^ ; ' paid to their proving method. Very recently, Ge and Jiang (GJ) reinvestigated previous derivation from a rigorously 
mathematic point of view They pointed out that HS misused the FK formula since it is usually applied for 
Ch ' Kolmogorov backward equation rather than Kolmogorov forward (or Fokker- Planck) equation. In fact, Chetrite 
I , and Gawedzki (CG) [ll| have given correct proof of the JE using the FK formula slightly earlier. Interestingly, however, 
' these two works seem to be apparently distinct though they employed the same mathematic formula: GJ constructed 
^ , a simple time-invariable integral (the proof of theorem 2.6 in Ref. Q), whereas CG was based on time-reversal concept. 

In addition to the JE, there are serval analogous equahties found in past several years p^. [l3l. [l^ [isl. [l6l. [itI . 
All of them have a type of 



(N ■ (exp[-^]) = 1, (1) 

> ■ 

' where ^ is a functional of the stochastic trajectory of a stochastic system, and angular brackets denote an average over 
Cf^ , the ensemble of the trajectories. For instance, A is total entropy change along a trajectory [l^. These equalities were 

■ caUcd integral fluctuation theorem (IFTs) to distinguish the detailed fluctuation theorem (DFTs) [Tsl. [l9l. [20l. [2ll. [2^ . 
I Due to the similarity between the IFTs with the JE in formality, one may naturally consider whether these equalities 

■ or even more general versions could be derived by GJ's approach as well. On the other hand, it should be intriguing 
' to clarify the relationship between their approach and the others, e.g. time-reversal, which was not performed in 

ON the previous work. In this work, we p resent our results about these two questions. With the help of the FK and 
Cameron-Martin-Girsanov (CMC) |23l. [23| formulas, we obtain a generalized integral fluctuation theorem. We find it 
could be explained from the concept of time-reversal. 



II. A TIME-INVARIABLE INTEGRAL 

We consider a general A^-dimension stochastic system x = {x;}, z = 1, • • • ,N described by a stochastic differential 
equation (SDE) ^ 

dx{t) = A(x, t)dt + B3 (x, t)dW{t) (2) 

where dW is an A^-dimensional Wiener process, A = {Ai} denotes a A^-dimensional drift vector and, B^/^ is the 
square root of a NxN positive definite diffusion matrix B = {-By }. We use Ito's convention for stochastic differentials. 
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Rather than directly solving Eq. ([2|). one usually converts the SDE into two equivalent PDEs of the conditional 
probability p(x, i|x',i') (i > t'): the forward Fokker-Planck equation (FPE) dtp = C{x,t)p with 



and the backward FPE dfp = — >C+(x',t')p with 



£+(x,i) = A,(x,i)0,, + iBy(x,t)a,.a,^.. (4) 

The initial conditions in both cases are p(x, t|x',<) = (5(x — x'). The two operators C and are adjoint each other. 
Here we use Einstein's summation convention throughout this work unless explicitly stated. 

Based on the symbols and definitions above, we find that, if u(x, t) satisfies a partial differential equation 

dt'u{x,t') = -C+{x,t')u{x,t') (5) 
- /-i(x,t') [drfi^,t') - /:(x,t')/(x,<')]"(x,t') 
+ /-i(x,t') [/:„(x,t')5(x,t')-.9(x,i'K^(x,t')] w(x,t'), 

where /(x, t') and g{x,t') are arbitrary normalized smooth positive functions [1^, Ca and are adjoint operators 
like Eqs. ^ and ([?]), which may be the same or different with those of the system we focus on, we have 

^y'dx/(x,i'Mx,t')=0. (6) 

The proof is obvious if one makes use of the adjoint property of these operators and notes that the derivative of the 
integral with respect to time is 



- J dx/(x, t')u{x, t') ^ Jdx [fdM^, t') + udrf{x, t')] . (7) 



This is simply generalization of GJ's idea, where the last three terms in Eq. ^ were absent [23] . Multiplying both 
sides of the PDF by /(x, i), one may see that this equation has a certain symmetry with respect to the functions / 
and u. In the following, we investigate Eq. (O with Liouvill- and Fokker-Planck- type La, respectively. Although one 
may think that the former is only a specific case of the latter (vanishing diffusion matrix) , the Liouvill-type may be 
more intriguing that will be seen shortly. 



III. LIOUVILL-TYPE Ca{y^,t) 



We assume 



£,g(x,i') = 29,,5,(x,t'), 



(8) 



and S = {5'i(x, t)} is an arbitrary vector having natural boundary condition. The coefficient 2 is for convenience in 
discussion. We rewrite Eq. ([5]) 

5t'u(x,t') = -/:+(x,t')w(x,t') (9) 

- /-i(x, t') [dv.f{x, t') - /:(x, t')/(x, t')] u{x, t') 

+ 2/-i(x, t') [a,.5.(x, t') + 5,(x, i')9.J u{x, t'). 
Using the FK and CMG formulas [l^l, we obtain a stochastic representation of the solution u{x,t) of the PDF given 

by 



u(x,t') = E^^'' {e-'^[-^'«'''(-)lg[x(t)]} 
where (j(x) is the final condition of u(x, t), and the functional 



(10) 



dr 



(11) 
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where the expectation E^'* is an average over all trajectories x(-) determined by Eq. ([2]) taken conditioned on 
x(i') = X, and the last term is Ito stochastic integral. Combining Eqs. (O and pU]) . we have 

J dx7(x', 0)£:^''° |e-'^W'^'''( )lg[x(i)] | = ^ dxf{x, i)g(x), (12) 

by choosing <' = 0. This equation has the same form as Eq. ^ for q{x.) = 1. We call it generalized integral fluctuation 
theorem (GIFT) in this work. Note other functions of q{x.) is still useful in practice. 



A. Relationship between GIFT and existing IFTs 

Existing several IFTs are specific cases of Eq. p^ . First case is to choose Si = and / = p^^{x, t), where p*^(x, t') 
is transient steady-state solution 



Then the functional (|lip is simply 



/:(x,i')p''(x,t') = 



J= - f 9rlnp""[x(T),T]dT. 
Jo 



(13) 



(14) 



GJ has analyzed this case with a specific final condition g(x) = 1 in great details. According to the system driven 
by a time-dependent conservative (a gradient of a potential) or non-conservative force, Eq. (jl2p reduces to Jarzynski 
equality 0, [![ and Hatano-Sasa relation [3| , respectively [§] . Here we are not ready to repeat the same derivatives 
but only to point out the necessary of the general final condition: one obtains the key Eq. (4) in the HS's work only 
choosing g(x) = 5(x — z). 

The second case is to choose Si to be the probability current Ji[f] of arbitrary normalized function /(x, t), where 



J,[/(x,t)] = A,(x,t)/(x,t) - -a,^.i3,,/(x,t). 
Substituting the current into Eq. (jlip and doing some derivations (details see the Appendix), we obtain 

/[x'(t'),t'] 



J = ln- 



./[x(t),t] 



2(S)^^ A; {B-^).. [xj - i(B-^)k.5.,(B-^)T 



dr 



(15) 



(16) 



where Ai = Ai — dxiBii/2 and xj = dxj/dT. Note that we used the Stratonovich' integral (denoted by "S" therein). 
Particularly, if the diffusion matrix is constant that is usually assumed in Langevin equation and /(x, t) is the solution 
of the stochastic system with a initial condition /(x',0), the functional is just the total entropy change Astot: the 
first term in Eq. (|16p is the system entropy change, and the second is the entropy change of environment p^ . Hence, 
Eq. inn is the IFT of the total entropy [ll, in this case. 

The last case is for the system being a nonequilibrium steady-state with a distribution p^^{x.Y_We choose / — p^^{x) 
and Si = Ji{p^^). Considering that d^iSi = 0, the functional becomes the housekeeping heat [l3| 



i(B-^),,a.,(B-^);^ 



dT 



(17) 



Of course, we can obtain this functional from Eq. (|16[) directly. Importantly, this IFT is even correct for time-dependent 
case [l3l. Although Eq. (fT^ does not include this equality, one can indeed prove it using CMC formula [ll[. 



IV. FOKKER-PLANCK-TYPE £a(x,t') 



If the diffusion matrix of La is not zero, we can still obtain a GIFT. Given 



/:a(x,<') 



-a,,a*(x,t') + i5,,a,>„(x,t') 



(18) 



Assuming a' = a5^(x, t) — 2a2(x, <:), we split this operator into a sum of a Fokker- Planck- and Liouvill operators, i.e. 

£a(x,i')= [>C;i(x,t') + 29,.4], (19) 



4 



where C\ is the same with Eq. ([18]) except that a' is replaced by a\. Substituting Eq. ((T9)) into Eq. ([5]) and doing a 
simple rearrangement, one may obtain Eq. ^ again if we redefine the drift vector and diffusion matrix of the operator 
C{yi,t') to be 

A', = ^,(x,t') + ai(x,<').9(x,t')//(x,t') 

Bl^ = i3,,(x,f')+6,,(x,t').9(x,t')//(x,t'), (20) 

and Si = a2(x, t')(7(x, t'). Hence, we have a GIFT (|12|) except that the expectation i?^^* is now an average over the 
trajectories generated from a new stochastic system 

dx(t) = A'(x, t)dt + B'3 (x, t)dW{t). (21) 

From this aspect, Eq. ([6]) with the Fokkcr-Planck-type Ca seems to not provide additional information compared to 
the Liouvill-type [28l |. 



V. GIFT AND TIME REVERSAL 



We must emphasize that the GIFT (fT2|) is derived mathematically (adjoint property of the FK operators) without 
resorting to any physical reason. It would be intriguing to establish a relationship between the equality and time- 
reversal like that done by Crook [l^ in interpreting JE [1, Q . CG recently suggested that the different fluctuation 
relations in literature may be traced to different time reversals [llj. In the remainder of the work, we use the same 
viewpoint to understand the origin of Eq. ([5]) and the GIFT. 

We first give a definition of a time reversal [ll[ . The variable Xi of the stochastic system is even or odd according 
to their rules under time reversal t' ^ t — t': if Xi yi = +Xi is even and Xi ^ yi = —xt is odd. We write 
them in abbreviation Xi ^ yi = EiXi and Ei = ±1. The drift vector splits into "irreversible" and "reversible" parts, 
A = A"'"' + A™^. Under a time reversal, these vectors are transformed into A = A'" + A''°'^, where 

ir(x,<') = e,A^^(y,s), (22) 
ir(x,t') = -e.,Ar{y,s), (23) 

y = {Vi} and s = t — t' {0 < t' < t). We also define a transformation of the diffusion matrix to be 

B,,(x,t') = e,e,Sy(y,s). (24) 

Note that no summation over repeated indices here. Because the splitting is completely arbitrary, one may think it 
to be an equivalent definition of a certain time reversal. 

Considering a forward FPE dspiy, s) = C{y, s)p(y, s) with the time-reversed drift vector and diffusion matrix. 



'C(y, s) = -dy,A,{y,s) + -dy,dy^Bij{y,s). (25) 



Substituting 

p(y,s)==/(x,i')t^(x,t') (26) 
into above equation and doing some simple derivations, we obtain 

at-«(x,i') = -/:+(x,t')i'(x,i') (27) 

- /-i(x, t') [9t,/(x, t') - £(x, t')/(x, t')] z;(x, t') 

+ 2/-i(x, t') [d^^sru) + sr{f)d^:\ v{^, f), 

where an irreversible probability current is defined to be 

STU) = ArMfM - ia,,B,fe(x,i')/(x,t') (28) 



We immediately see that Eq. © would be the same with Eq. (P7|) if Si = Sl"{f). Under this circumstance, Eq. 

is a obvious consequence of the specific time reversal of the stochastic system in that the integral of the left-hand side 
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of Eq. ([6]) equals to J dyp{y, s), which is conservative with respect to time. On the other hand, considering that the 
sphtting of the drift vector is completely arbitrary, we can define a time reversal for a given vector Si by 



/(x,t') 



A" 



c,t') = A,(x,t')-A"-(x,t') 



(29) 
(30) 



Then the GIFT for any vector S is a consequence of the new defined time reversal of the stochastic system. The 
reader should be reminded that these time reversals are not always meaningful or easily realizable in physics. Before 
ending the section, we give two comments about Eq. (j26p . First, for a homogenous diffusion process that satisfies 
the detailed balance, this equation has been used earlier to connect the forward FK solution J3(x, t) and backward 
FK solution w(x, t), where f{x,t) 



[Eq. pU] are crucial for achieving DFTs [1 
therein) [ll|- 




Second, Eqs. ([26|) and the stochastic representation solution w(x, t') 
. CG has given a detailed discussion about this point (Proposition 1 



VI. CONCLUSION 

In this work, we present a generalized integral fluctuation theorem for diffusion processes using the famous FK 
and CMG formulas. Although one might think that this result (|12|) is almost trivial from the point of view of 
time-reversal that we explained here, the derivation by constructing a time-invariable integral is novel. We close this 
work by pointing out several directions that we can further improve the current results. First is to give rigorous 
conditions for applying FK and CMG formulas. Second, we may extend the results for the continuous diffusion 
processes to discrete state case, e.g. chemical master equation. Finally, the most interesting to us is to apply the 
GIFT in practical problems, e.g. free energy reconstruction of single-molecule manipulation. 

We appreciate Dr. Chetrite for sending their intriguing work to us. The discussion about the connection between 
the GIFT and time-reversal is mainly inspired by them. This work was supported in part by Tsinghua Basic Research 
Foundation and by the National Science Foundation of China under Grant No. 10704045 and No. 10547002. 



APPENDIX: DERIVATION OF THE TOTAL ENTROPY 

Substituting Eq. (|15p into Eq. (jlip and separating the terms into /-dependent and -independent parts, we obtain 

/•* r 11 

dT + {d^Anf){B^),jdWj (31) 



drlnf + AAi In / + -B.jd^^d^. In / 



Employing the SDE the second integral in above equation becomes 



5r In / + -Bijdxidx^ In / 



dr + dx, \nfdx,, 



(32) 



which is just the system entropy change J^, din / = In [f {x{t) , t) / f {x{t') , t')] according to the Ito formula The 
entropy change of environment in Eq. (jl6p is derived from the first integral in Eq. (j3ip by converting the Ito's stochastic 
integral into the Stratonovich's: 



2MB-^).,,dW, = (S) / 2Ai(B-5)ydWj - l{Bi)A 



2(B-^)jAi 



dr. 



(33) 



An alternative method is based on time-reversal technique. Previous discussion shows that for u(x, t') satisfying 
Eq. ^ and Si ~ Ji{f), we can define a time-reversal with A"'' ~ A and A™^ = 0. Then f{x,t')u{x,t') = p{y,s) 
is a solution of the FPE p5|) . Because / is arbitrary, we may choose / = 1 and w{x,t') = p{y,s). Hence we have 
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w(x, t') = ?«(x, t') exp[— ln/(x, i')]. Obviously, the functional of the stochastic representation of ^(x, t') is the first 
integral of Eq. ([?T|) . Like the proof of the GIFT, this derivation shows the power of time- reversal technique again. 
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